Domain formation in transitions with noise and time-dependent bifurcation parameter 
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The characteristic size for spatial structure, that emerges when the bifurcation parameter in model partial 
differential equations is slowly increased through its critical value, depends logarithmically on the size of added 
noise. Numerics and analysis are presented for the real Ginzburg-Landau and Swift-Hohenberg equations. 
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Many physical systems undergo a transition from a 
spatially uniform state to one of lower symmetry. Clas- 
sical examples are the formation of magnetic domains 
and the Rayleigh-Benard instability [1]. Such systems 
are commonly modeled by a simple differential equation, 
having a bifurcation parameter with a critical value at 
which the spatially uniform state loses stability. Noise is 
often assumed to provide the initial symmetry-breaking 
perturbation permitting the system to choose one of the 
available lower-symmetry states, but is not often explic- 
itly included in mathematical models. However, when 
the bifurcation parameter is slowly increased through its 
critical value it is necessary to consider noise explicitly. 

The phenomenon of delayed bifurcation and its sen- 
sitivity to noise has been reported in the case of non- 
autonomous stochastic ordinary differential equations [2] ; 
here the corresponding phenomenon is examined in par- 
tial differential equations. A characteristic length for the 
spatial pattern is demonstrated from a stochastic par- 
tial differential equation (SPDE), supported by numeri- 
cal simulations. Noise is added in such a way that it has 
no correlation length of its own (white in space and time) 
and a finite difference algorithm is used whose continuum 
limit is an SPDE. 

The mathematical description of transitions is in terms 
of a space-dependent order parameter Y and a bifurca- 
tion parameter g. Because it is the simplest model with 
the essential features, the real Ginzburg-Landau equa- 
tion (GL) is considered first. Results are also presented 
for the Swift-Hohenberg equation (SH), that is more ex- 
plicitly designed to model Rayleigh-Benard convection. 

When the bifurcation parameter g is constant the fol- 
lowing is found. For g < 0, in both GL and SH, the 
solution with Y everywhere is stable. In GL for g > 
one sees a pattern of regions where Y is positive and 
regions where Y is negative (domains) separated by nar- 
row transition layers. In SH for g > there is a structure 
resembling a pattern of parallel rolls, interrupted by de- 
fects. 

When g is slowly increased through in the presence of 
noise a characteristic length is produced as follows. The 
field Y remains everywhere small until well after g passes 
through 0. At g ~ g c , where 



\i is the rate of increase of g and e is the amplitude of 
the noise, Y at last becomes 0(1) and the spatial pattern 
present is frozen in by the nonlinearity. Thereafter one 
observes spatial structure with characteristic size propor- 
tional to (| lne|//i)4. In GL this length is the typical size 
of the domains; in SH it is the typical distance betwen 
defects. 

The results reported here were obtained by solv- 
ing SPDEs [3] of the following dimensionless form for 
stochastic processes Y depending on x and t: 



dY = [g(t)Y - Y 3 + £Y]dt + edW. 



(2) 



(1) 



The equations were solved as initial value problems, 
with g(t) = jit slowly increased from —1 to 1. Here 
Y : [0,L] m x ±] x n -> K, Q is a probability space 
and W is the Brownian sheet [4] , the generalisation of the 
Wiener process (standard Brownian motion) to processes 
dependent on both space and time. Periodic boundaries 
in x are used so that any spatial structure is not a bound- 
ary effect. The constants /i, e and are all <C 1. Results 
are reported for C = A (GL) and C = -(1 + A) 2 (SH) 
where A = YliLi the Laplacian in lZ m . 

In the first order finite difference algorithm for numer- 
ical realisations of the lattice version of (2), yt+At{i) is 
generated from y t (i) as follows: 

Vt+At(i) = Vt(i) + \}ityt{i) - y 3 (i) + Cy t (i)]At 

-he (Ax)~^ n t (i)VAt. (3) 

In (3), yt(i) is numerical approximation to the value of Y 
at site i at time t and C is the discrete version of C. The 
rtt(i) are Gaussian random variables with unit variance, 
independent of each other, of the values at other sites, 
and of the values at other times. 

It is also possible to introduce multiplicative noise, for 
example to make g a random function of space and time 
[5,6]. The effect in that case is proportional to the mag- 
nitude of the noise and is thus less dramatic at small 
intensities than that of additive noise. 

The timing of the emergence of spatial structure can be 
understood by deriving the stochastic ordinary differen- 
tial equation for the most unstable Fourier mode, which 
is of the form [7] 
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dy= [g(t)y-y 3 ]dt + edw, 



(4) 



where w is the Wiener process. Trajectories of (4) re- 
main close to y = until well after g = 0, then jump 
abruptly towards one of the new attractors (Figure 1). 
The value of g at the jump can be determined by solving 
the linearised version ; for /i <C 1 it is a random variable 
with mean approximately g c and standard deviation pro- 
portional tO fjL [8] . 
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FIG. 1. Dynamic pitchfork bifurcation with noise. The 
dotted lines are the loci of stable fixed points of y = gy — y 3 
as a function of g. The solid lines are solutions of the 
non- autonomous SDE (4) with g — fit, for noise levels 
e = 1(T 3 , 1(T 6 , 1(T 9 , 1(T 12 . (In each case fi = 0.01 and 
the initial condition is y = 1.0 at g = —1.0.) 



The Ginzburg-Landau equation is a simple model of 
a spatially extended system where a uniform state loses 
stability to a collection of non-symmetric states. When 
g is fixed and positive in this equation, a pattern of do- 
mains is found. In each domain, Y is close either to 
y/g or to — yfg. The gradual merging of domains on ex- 
tremely long timescales [9] is not the subject of this pa- 
per; here the focus is on how the domains are formed by 
a slow increase of the bifurcation parameter through 0. 
An example is depicted in Figure 2: a pattern of domains 
emerges when Y is everywhere small and is frozen in at 
g ~ g c . When Y is small an excellent approximation to 
the correlation function, c(x) = (Y t (v)Y t (x + 1>)}, can be 
calculated from the solution of the linearised version of 
(2) (that is, without the cubic term). The correlation 
length at g = g c becomes the characteristic length for 
spatial structure after g = g c . 

For GL, the solution of the linearised version of (2) is: 

Y t (x)= [ G(t,-^,v)f(v)dv + 

J[0,L] m 

e [ [ G(t,8,x,v)dvdW a (v), (5) 

J-l J[0,L]™ 
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FIG. 2. Dynamic transition, GL, one space dimen- 
sion. Four configurations, Yt(x), are shown from one nu- 
merically-generated realisation of the SPDE (note the differ- 
ent vertical scales). Nonlinear terms become important when 
g ~ 0.64; their effect is to freeze in the spatial structure. 
(L = 300, /x = 0.01, e= 10" 10 .) 



where G(t,s,x,v) 



[4tt(£-s)] ™ exp (- 



4(t-s) 



2 ) 



with x — v understood modulo [0,L] m . The first term, 
dependent on the initial data /(#), relaxes quickly to 
very small values and remains negligible if 2fi\ lne| < 1. 
The correlation function is therefore obtained from the 
second integral in (5). The mean of the product of two 
such stochastic integrals is an ordinary integral [4] . Per- 
forming the integration over space [7], assuming that 

L > (f )*> § ives 
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c(x) = e 5 



t e /i(t 2 -S 2 ) e8(t _ s) 



-ds. 



(6) 



Before g approaches 0, the correlation function differs 
by only O(j^) from its static (g =constant) form [7]; 
it remains well-behaved as g passes through and, for 
g > y/ji, is we U approximated by: 



c(x) 



eV 4 



(7) 



For -j= < g < g c , typical values ofY(x) increase expo- 
nentially fast and the correlation length is proportional 
to y/t. Effectively noise acts for g < y/Jl to provide an 
initial condition for the subsequent evolution. At a value 
of g that is a random variable with mean g ~ g c and 
standard deviation proportional to /i, the cubic nonlin- 
earity becomes important. Its effect is to freeze in the 
spatial structure; no perceptible changes occur between 
g = g c and g = 1. 

In one space dimension it is possible to put the scenario 
just described to quantitative test by producing numer- 
ous realisations like that of Fig. 2 and recording the num- 
ber of times Y crosses upwards through in the domain 
[0, L] at g = 1. In Fig. 3 the average number of upcross- 
ings is displayed as a function of the sweep rate /i. The 
solid line is the expected number of upcrossings of zero, 
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for a field with correlation function (7) at g = g c [10]. 
The hypothesis that the spatial pattern does not change 
after g = g c is succesful. 
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FIG. 3. Number of zero crossings after a dynamic transi- 
tion. The dots are the mean number of upcrossings of at 
g = 1 in numerical realisations of GL in one space dimen- 
sion. The solid line is the prediction based on the assumption 
that the correlation function (7) is valid until g = g c , after 
which time the spatial pattern does not change, (e = 10 ~ 4 
and L = 800.) 



In one space dimension, the solution of the SPDE (2) 
is a stochastic process with values in a space of contin- 
uous functions [3,12,13]. That is, for fixed uo G and 
t G [— ^,^], one obtains a configuration, Y t (x), that is 
a continuous function of x. This can be pictured as the 
shape of a string at time t that is constantly subject to 
small random impulses all along its length. In more than 
one space dimension, however, the Y t (x) are not neces- 
sarily continuous functions but only distributions [3,12]. 
Typically the correlation function c(x) diverges at x = 0. 
In the dynamic equations studied here, however, the di- 
vergent part does not grow exponentially for g > 0, and 
by 9 = 9c it is only apparent on extremely small scales, 
beyond the resolution of any feasible finite difference al- 
gorithm. Figure 4 depicts configurations at g = 1 from 
realisations of (2) in two space dimensions. In Figures 
4(a) and 4(b) (GL) one sees that a faster rate of increase 
of g results in a smaller average domain size. The SPDEs 
were simulated on a grid of512x512 points with second 
order timestepping [13]. 




FIG. 4. Two-dimensional pattern at g — 1: smaller u 
means larger characteristic length. In black regions Y < 0; in 
white or grey regions Y > 0. In GL, (a) and (b), the typical 
domain size decreases with u , the rate of increase of g. In SH, 
(c) and (d), where there is a short-range structure resembling 
parallel rolls, the effect of reducing u is to reduce the number 
of defects. 4(a): GL, L = 300, e = 10" 5 , a = 0.03. 4(b): GL, 
L = 300, e = 10~ 5 , a = 0.003. 4(c): SH, L = 200, e = 10~ 5 , 
u - 0.01. 4(d): SH, L = 200, e = 10~ 5 , u = 0.001. 

The essential difference between the Swift-Hohenberg 
and Ginzburg-Landau models is that the first spatial 
Fourier mode to become unstable has k = 1 rather than 
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k = 0. Hence there is a preferred small-scale pattern 
that resembles the parallel rolls seen in experiments. 
However, there is no preferred orientation of the roll pat- 
tern and when the correlation length is smaller than the 
system size, many defects are found, separating regions 
where the rolls have different orientations. When g is 
increased through 0, the number of defects resulting de- 
creases when fi decreases - Fig. 4(c) and (d). Here a grid 
of 300 x 300 points was used with first order timestep- 
ping. 

A notable feature of dynamic bifurcations and dynamic 
transitions is that the evolution for g > is independent 
of the initial conditions (provided they are such that 
that the system descends into the noise). Noise acts, 
near g = 0, 



to wipe out the memory of the system and to provide 
an initial condition for the subsequent evolution. The 
correlation function (7) is, for example, a natural initial 
condition for studying the dynamics of defects and phase 
separation because it emerges from a slow increase to su- 
percritical of the bifurcation parameter in the presence 
of space-time noise, mimicking an idealised experimen- 
tal situation. 

In summary, dynamic transitions are analysed in models 
of spatially extended systems with white noise. The 
correlation length that emerges from the noise during a 
slow sweep past g = is frozen in by the nonlinearity as 
a characteristic length proportional to (| lne|//x)4 where 
H is the rate of increase of the bifurcation parameter and 
e is the amplitude of the noise. 
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